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[i]  We  present  the  first  SHIMMER  observations  of  the 
diurnal  variation  of  mesospheric  hydroxyl  (OH).  We 
compare  our  data  with  Aura  Microwave  Limb  Sounder 
(MLS)  observations  at  about  13h  local  time  near  55°N  and 
find  very  good  agreement.  This  validates  the  Spatial 
Heterodyne  Spectroscopy  technique  for  space-borne 
optical  remote  sensing  applications.  We  extend  our 
analysis  to  other  local  times,  not  observed  by  MLS,  for 
latitudes  near  55°N  in  the  summer  of  2007.  At  74  km,  we 
find  excellent  agreement  with  a  photochemical  model,  but 
above  76  km,  significant  model/data  differences  in  the 
shape  of  the  OH  diurnal  variation  are  observed. 
Citation:  Englert,  C.  R.,  M.  H.  Stevens,  D.  E.  Siskind,  J.  M. 
Harlander,  F.  L.  Roesler,  H.  M.  Pickett,  C.  von  Savigny,  and  A.  J. 
Kochenash  (2008),  First  results  from  the  Spatial  Heterodyne 
Imager  for  Mesospheric  Radicals  (SHIMMER):  Diurnal  variation 
of  mesospheric  hydroxyl,  Geophys.  Res.  Lett.,  35,  L19813, 
doi:  1 0. 1 029/2008GL035420. 

1.  Introduction 

[2]  In  the  mesosphere,  hydroxyl  (OH)  is  predominantly 
produced  by  the  photodissociation  of  water  vapor  (H20) 
and  its  use  as  a  proxy  for  mesospheric  H20  has  been 
demonstrated  [Summers  et  al.,  1997,  2001],  Measurements 
of  mesospheric  hydroxyl  (OH)  have  been  performed  from 
aircraft,  balloons,  and  satellites.  The  aircraft  measurements, 
using  far  infrared  emission,  were  limited  to  a  vertical 
resolution  of  about  25  km  in  the  mesosphere  [Englert  et 
al.,  2000].  Balloon  measurements  have  provided  column 
densities  above  about  40  km  [Pickett  et  al.,  2006a],  Recent¬ 
ly,  only  four  instruments  in  Earth  orbit  have  been  capable  of 
providing  high  vertical  resolution  OH  altitude  profiles  in  the 
mesosphere:  The  Imaging  Spectrometric  Observatory  (ISO) 
[Morgan  et  al.,  1993],  the  Middle  Atmospheric  High 
Resolution  Spectrograph  Investigation  (MAHRSI)  [Conway 
et  al.,  1999],  the  Optical  Spectrograph  and  Infrared  Imager 
System  (OSIRIS)  [Gattinger  et  al.,  2006],  and  the  Micro¬ 
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wave  Limb  Sounder  (MLS)  on  Aura  [Pickett  et  al. ,  2006a, 
2006b,  2008], 

[3]  ISO  OH  data  was  reported  for  one  day,  one  latitude, 
and  one  local  time.  The  MAHRSI  data  are  limited  to  two 
Space  Shuttle  flights  in  1994  and  1997  which  provided  very 
little  local  time  (LT)  coverage  for  any  given  latitude  and 
both  OSIRIS  and  MLS  observations  are  made  from  sun 
synchronous  platforms  that  essentially  sample  only  two 
local  times.  In  contrast,  the  Spatial  Heterodyne  Imager  for 
Mesospheric  Radicals  (SHIMMER)  measures  mesospheric 
OH  over  the  diurnal  cycle  and  the  first  results  are  discussed 
here.  Like  the  diurnal  measurements  of  H02,  O3  and  H20, 
SHIMMER  data  are  expected  to  be  an  invaluable  tool  to 
assess  our  understanding  of  the  processes  describing  meso¬ 
spheric  HOx  (H02,  OH,  H)  [Sandor  and  Clancy,  1998]. 

[4]  SHIMMER  was  launched  on  board  the  Space  Test 
Program  Satellite-1  (STPSat-1)  on  March  9,  2007.  It  meas¬ 
ures  the  OH  solar  resonance  fluorescence  by  observing  the 
A2£+-X2n(0,0)  band  near  309  rnn.  SHIMMER  uses  the 
innovative  Spatial  Heterodyne  Spectroscopy  technique 
which  offers  considerable  advantages  in  reducing  size  and 
weight  compared  to  conventional  grating  spectrographs 
[Harlander  et  al.,  2003],  SHIMMER  has  no  moving  optical 
components  and  images  the  Earth’s  limb,  rather  than  scan¬ 
ning  it.  The  STPSat-1  orbital  inclination  and  orbit  altitude 
result  in  observations  up  to  58°  latitude  with  a  local  time 
precession  of  ~30  minutes  per  day,  allowing  SHIMMER  to 
measure  over  the  diurnal  cycle  [Stevens  et  al.,  2008], 

[5]  In  this  study,  we  first  compare  the  SHIMMER  OH 
measurements  with  MLS  data  at  the  local  time  of  the  MLS 
measurements  (~13h).  Subsequently,  we  investigate  the 
local  time  variation  of  OH  as  measured  by  SHIMMER  in 
the  upper  mesosphere  and  compare  to  results  from  a  one¬ 
dimensional  photochemical  model. 

2.  SHIMMER  Observations 

[6]  The  STPSat-1  satellite  is  in  a  35.4°  orbital  inclination 
and  the  SHIMMER  field  of  view  is  oriented  perpendicular 
to  the  satellite  velocity  direction  towards  the  summer  pole 
so  that  it  observes  to  about  58°  latitude.  It  typically  images 
the  limb  between  ~33-97  km  altitude,  with  2  km  altitude 
sampling,  a  12  second  integration  time,  and  a  sampling 
period  of  20  seconds.  SHIMMER  measurements  provide 
very  dense  coverage  within  the  53°N-58°N  latitude  band  in 
the  northern  summer,  with  over  500  limb  images  per  day 
[Stevens  et  al.,  2008].  We  therefore  focus  on  this  latitude 
region  for  our  comparison  with  MLS  data  and  model 
results. 

[7]  The  raw  SHIMMER  interferograms  are  downlinked 
and  processed  on  the  ground.  The  dark  field  correction  is 
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Figure  1.  (top)  Averaged  SHIMMER  spectrum  at  72  km 
tangent  altitude  (black).  Fitted  spectrum  including  the 
Rayleigh  scattered  solar  background  and  the  OH  emission 
lines  (red).  Estimated  Rayleigh  scattered  background  (blue), 
(bottom)  Difference  between  the  measured  spectrum  and  the 
Rayleigh  scattered  background  (black).  Fitted  OH  emission 
lines  (red). 

performed  using  on-orbit  measurements,  the  flat-field  cor¬ 
rection  is  performed  applying  the  unbalanced  arm  flat- 
fielding  approach  [Englert  and  Harlander,  2006]  based  on 
laboratory  measurements,  and  the  phase  correction  assumes 
a  low  resolution  phase  shift  [ Englert  et  al.,  2004].  Further¬ 
more,  corrections  for  cosmic  ray  effects  on  the  detector, 
thermal  defocus,  and  grating  imperfections  are  applied.  The 
radiometric  calibration  is  based  on  laboratory  measurements 
using  a  National  Institute  of  Standards  and  Technology 
(NIST)  traceable  UV  light  source.  A  representative,  aver¬ 
aged  spectrum  for  a  tangent  point  height  of  72  km  is  shown 
as  the  black  line  in  Figure  1  (top). 

[8]  The  separation  of  the  OH  resonance  fluorescence 
signal  from  the  Rayleigh  scattered  background  is  performed 
using  a  method  that  is  similar  to  the  one  used  for  the 
MAHRSI  data  [Conway  et  al.,  1999].  This  includes  fitting 
the  spectral  shapes  of  the  OH  emission,  the  Rayleigh 
scattered  solar  flux,  a  baseline  offset,  and  the  instrumental 
line  shape  function.  The  solar  background  shape  was 
obtained  from  on-orbit  SHIMMER  measurements  of  direct 
sunlight.  A  limb  observation  at  72  km,  the  estimated 
background  and  the  best  fit  of  the  total  signal  are  shown  in 
Figure  1  (top).  The  measured  signal  after  the  subtraction  of 
the  background  and  the  fitted  OH  emission  spectrum,  con¬ 
volved  with  the  instrumental  line  shape  function  are  shown  in 
the  bottom  panel.  We  determine  the  brightness  of  the  OH 
emission  for  each  altitude  using  this  proven  approach. 

[9]  The  altitude  registration  of  the  limb  images  is 
achieved  using  the  satellite  position  and  attitude  data  as 
measured  by  the  on-board  Global  Positioning  System  and 
star  tracker.  SHIMMER  does  not  have  the  capability  to 
verily  the  pointing  accuracy  with  on-orbit  observations  of 
known  targets  like  stars.  Therefore,  in  order  to  correct  for 
any  systematic  errors  in  the  alignment  of  the  star  tracker 
with  respect  to  SHIMMER,  including  the  on-orbit  thermo¬ 
elastic  distortion  of  the  satellite,  we  employ  the  Tangent 
Height  Retrieval  by  UV-B  Exploitation  (TRUE)  method, 


which  has  previously  demonstrated  an  accuracy  of  better 
than  1  km  [von  Savigny  et  al.,  2005].  Using  a  subset  of 
SHIMMER  data  used  for  this  work,  we  find  an  altitude 
offset  of  2.81  ±  0.35  km,  which  we  apply  to  the  data. 

[10]  To  retrieve  the  OH  altitude  concentration  profile 
from  the  OH  brightness  profile  we  use  a  Twomey  regular¬ 
ized  least  squares  method  as  described  by  Conway  et  al. 
[1999].  Our  version  is  nearly  identical  to  the  algorithm  that 
was  used  for  the  MAHRSI  data  analysis  [Summers  et  al., 
1997,  2001;  Conway  et  al.,  1999,  2000],  including  the  same 
OH(0,0)  emission  rate  factors  (g  factors)  [Stevens  and 
Conway,  1999].  The  main  difference  is  that  the  analysis  is 
adjusted  for  the  slightly  smaller  SHIMMER  passband. 

3.  Comparison  of  SHIMMER  and  MLS  OH 
Observations 

[11]  To  compare  SHIMMER  OH  results  with  another 
coincident  measurement,  we  use  observations  from  the 
MLS  instrument  on  NASA’s  Aura  satellite.  We  choose  data 
from  July  16,  2007  during  which  the  two  instruments  were 
observing  the  limb  at  the  same  local  time.  We  take  a  zonal 
average  of  SHIMMER  data  (version  1.0)  at  12.6  ±  2.6  hours 
LT  and  a  zonal  average  of  MLS  data  (version  2.21)  at 
12.8  hours  LT  between  53-58°N.  Figure  2  illustrates  the 
comparison  of  the  observations  between  60-90  km  altitude 
with  shaded  envelopes  representing  the  combined  lcr  sta¬ 
tistical  and  systematic  uncertainties.  The  SHIMMER  uncer¬ 
tainty  is  the  root  sum  square  (rss)  of  an  altitude  dependent 
statistical  error  derived  from  the  baseline  variation  in  the 
spectrum,  systematic  errors  that  are  predominantly  due  to 
the  correction  of  grating  imperfections  and  a  7%  rss 
uncertainty  in  the  OH  g  factors  [Conway  et  al,  1999], 
The  MLS  uncertainty  includes  the  statistical  error  and  a 
10%  systematic  error.  The  agreement  is  excellent,  especially 
considering  that  the  two  observations  use  completely  dif¬ 
ferent  remote  sensing  techniques. 


OH  Density  (107  cm'3) 


Figure  2.  Comparison  of  OH  vertical  density  profiles 
measured  by  SHIMMER  (black)  and  MLS  (blue)  and 
calculated  by  the  NRL-CHEM1D  photochemical  model 
(red).  The  model  result  is  multiplied  with  the  SHIMMER 
weighting  functions  to  match  the  altitude  resolution.  The 
model  result  is  initialized  with  MLS  H20  for  the  same  day, 
latitudes  and  local  time. 
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Figure  3.  Water  vapor  profiles  measured  by  MLS  for  July 
16,  2007  and  June  30,  2007.  The  profiles  show  only  small 
differences  (less  than  ±10%)  between  60-85  km. 

[12]  Also  shown  in  Figure  2  is  the  result  of  an  NRL- 
CHEM1D  photochemical  model  calculation  [ Summers  et 
al.,  1997,  2001]  using  standard  JPL06  reaction  rates 
[Sander  et  al.,  2006],  MLS  FLO  (Version  2.2)  [Lambert  et 
al.,  2007],  temperature,  and  pressure  profiles  for  the  same 
latitude  range.  The  zonally  averaged  MLS  water  vapor 
profile  for  July  16  used  in  the  model  calculation  is  shown 
in  red  in  Figure  3.  SHIMMER  OH  data  between  60-90  km 
are  in  good  agreement  with  the  model.  This  is  in  contrast  to 
MAHRSI  observations  of  solar  resonance  fluorescence, 
which  indicated  a  25-30%  model  overprediction  in  the 
mesosphere  [Conway  et  al.,  1999,  2000]. 

[13]  It  is  important  to  note  that  while  the  averaging 
kernels  for  SHIMMER  OH  and  MLS  OH  are  3-7  km  wide 
between  65-90  km,  the  averaging  kernels  for  the  MLS  H20 
are  between  11-16  km  wide  [Lambert  et  al.,  2007]. 

[14]  Figure  2  also  serves  as  a  verification  that  the 
SHIMMER  OH  measurement  does  not  suffer  from  a  sig¬ 
nificant  scattered  light  problem.  In  addition  to  the  good 
agreement  between  SHIMMER  and  MLS  results,  the  re¬ 
trieved  OH  density  approaches  zero  above  85km,  which 
confirms  that  OH  signal  from  below  does  not  contaminate 
the  higher  altitudes. 

4.  Mesospheric  OH  Over  the  Diurnal  Cycle 

[15]  The  local  time  precession  of  the  SHIMMER  meas¬ 
urements  allows  us  to  investigate  the  diurnal  variation  of 
OH,  which  we  accomplish  by  comparing  zonally  averaged 
SHIMMER  data  with  results  of  the  NRL-CHEM1D  model. 
We  use  zonally  averaged  SHIMMER  OH  profiles  between 
55-58°N  and  between  27  June  and  2  August,  2007,  during 
which  the  average  local  time  decreased  monotonically  from 
20.6  to  5.1  hours.  The  results  are  shown  in  Figure  4  for  six 
altitudes  between  74  km  and  84  km.  The  SHIMMER  points 
are  3  day  averages,  resulting  in  one  data  point  every  ~1.5h 
LT.  OH  observations  after  20. 6h  and  before  5. lh  LT  are  not 
shown  due  to  the  lack  of  sunlight  required  for  the  fluores¬ 
cence  observations. 

[16]  For  the  model  results,  we  use  zonally  averaged  MLS 
H20  profiles  for  the  same  latitude  region  and  time  periods 


measured  at  about  13h  LT.  The  individual  water  profiles 
show  only  small  changes  during  this  period  as  indicated  by 
the  two  profiles  for  June  30  and  July  16  in  Figure  3. 

[17]  At  74  km  the  agreement  between  the  SHIMMER  OH 
data  and  NRL-CHEM1D  is  excellent.  Here,  the  shape  of  the 
diurnal  variation  is  dominated  by  the  variation  of  the  OH 
production  by  photodissociation  of  water,  which  is  a  func¬ 
tion  of  the  solar  zenith  angle.  Therefore,  the  OH  concen¬ 
tration  peaks  around  local  noon,  for  the  smallest  solar  zenith 
angle. 

[is]  Above  76  km,  on  the  other  hand,  the  measured  shape 
of  the  diurnal  variation  is  significantly  different  from  the 
model  results.  At  these  altitudes,  the  modeled  OH  peaks 
later  in  the  afternoon,  resulting  from  the  interaction  of  the 
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Figure  4.  Diurnal  behavior  of  OH  as  measured  by 
SHIMMER  (black  diamonds)  compared  with  CHEM1D 
calculations  (red  squares).  The  error  envelope  of  the 
SHIMMER  measurements  reflects  an  estimate  of  the 
combined  lcr  systematic  and  statistical  uncertainties. 
SHIMMER  local  times  vary  from  20.6  h  on  27  June  (Day 
178)  to  5.1  h  on  2  August  (Day  214).  The  CHEM1D  curves 
represent  a  composite  of  12  separate  calculations  for  the 
dates  indicated  by  the  top  axis.  Each  calculation  uses  the 
zonal  average  MLS  H20  profile  at  ~  1 3  h  LT  for  the  same 
latitude  region  (55-58°N).  Each  red  square  represents  the 
model  result  for  the  local  time  given  by  the  bottom  axis.  The 
model  results  are  smoothed  with  the  SHIMMER  altitude 
resolution. 
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HOx  species  with  ozone  [e.g.,  Prather,  1981;  Summers  et 
al.,  2001],  SHIMMER,  however,  shows  a  systematically 
different  shape  of  the  diurnal  variation  with  absolute  density 
differences  of  two  or  more,  which  is  well  beyond  the 
uncertainties. 

5.  Discussion  and  Conclusions 

[19]  We  have  compared  mesospheric  OH  density  profiles 
of  SHIMMER  and  MLS  for  the  same  day,  the  same  latitudes 
and  at  the  same  local  time.  The  results  are  in  excellent 
agreement,  proving  that  the  innovative  SHS  optical  tech¬ 
nique  is  valuable  for  long  duration  space-bome  UV  remote 
sensing. 

[20]  We  have  also  compared  SHIMMER  and  MLS  OH 
observations  near  13LT  with  results  from  a  one-dimensional 
photochemical  model  using  standard  JPL06  chemistry, 
where  the  model  is  initialized  with  MLS  H2O.  We  find  that 
both  SHIMMER  and  MLS  data  are  in  good  agreement  with 
the  model  results. 

[21]  Finally,  we  have  presented  the  first  satellite  obser¬ 
vations  of  mesospheric  OH  over  nearly  the  entire  diurnal 
cycle  and  here  we  find  discrepancies  with  the  model  above 
76  km.  In  an  effort  to  reconcile  these  differences  we 
considered  the  reaction  rate  uncertainty  of  the  dominant 
reaction  which  partitions  OH  and  H  (H±03^0H±02) 
[Mlynczak  et  al.,  1998]  and  which  is  of  particular  impor¬ 
tance  in  the  altitude  region  where  we  see  the  differences.  We 
scaled  this  rate  by  ±50%  which  is  about  twice  the  uncer¬ 
tainty  given  by  Sander  et  al.  [2006]  and  found  that  the  OH 
around  80  km  responded  significantly,  but  it  did  not 
improve  the  model/data  comparison.  Furthermore,  diumally 
varying  temperatures  due  to  tides  increasingly  become  more 
important  at  the  higher  altitudes.  Likewise,  because  H20 
has  a  strong  vertical  gradient  at  80  km,  it  too  should  vary 
with  tidal  upwelling  and  downwelling;  however,  the  mag¬ 
nitude  of  the  OH  model-data  discrepancy  is  too  large  to  be 
remedied  by  a  simple  scaling  of  the  H20  model  input.  Also, 
the  HOx  chemical  lifetime  is  increasing  rapidly  above 
~80  km  [Brasseur  and  Solomon,  1986]  suggesting  that 
transport  effects  may  play  a  role.  Our  model  calculations 
shown  here  did  not  explicitly  account  for  diffusive  trans¬ 
port;  however,  a  test  simulation  with  vertical  eddy  diffusion 
(Kzz)  confirmed  that  it  can  be  of  importance  at  these 
altitudes.  Finally,  since  these  SHIMMER  measurements 
are  made  at  the  fringe  region  of  noctilucent  clouds,  the 
redistribution  of  H20  through  ice  particle  growth  and 
sublimation  might  be  relevant  [Murray  and  Plane,  2005], 

[22]  A  detailed  parametric  study  of  all  these  factors  is 
beyond  the  scope  of  the  present  paper.  We  suggest  that  the 
high  quality  of  the  SHIMMER  data  with  its  substantial 
diurnal  and  seasonal  coverage  (as  of  this  writing,  SHIMMER 
observations  have  covered  more  than  one  year),  combined 
with  other  observations  such  as  that  of  O3  and  temperature 
from  SABER  on  TIMED  [ Huang  et  al.,  2008]  would  be 
invaluable  in  providing  new  insights  about  the  roles  of 
chemistry  and  dynamics  in  the  upper  mesosphere  over  the 
diurnal  cycle. 
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